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1 INTRODUCTION 

IT) 

The 1-20 keV spectra of Galactic Black Hole Candidates 
(*] (BHC) are characterized by independently variable soft 
ClJ^ckbody-like components and power laws. Recent work 
ion high quality spectra from EXOSAT and HEAO-1 (Done 
£~et al. 1992), Tenma (Kitamoto et al. 1990), Ginga (Ebi- 
sawa 1991; Tanaka 1992; Inoue 1993; Ueda, Ebisawa & Done 
1994) and BBXRT (Marshall et al. 1993) have revealed spec- 
tral features above 6 keV which are probably due to absorp- 
tion and emission by iron. The strongest feature is an ab- 
sorption edge occurring above 7 keV which appears smeared 
relative to a simple photoelectric edge (e.g. Ebisawa 1991). A 
remarkable aspect of such an edge is that the limits on the 
equivalent width of any associated emission line are much 
lower than expected on the basis that about one third of all 
photons absorbed should result in fluorescent line photons 
around 6.4 keV (Done et al. 1992; Tanaka 1992; Marshall et 
al. 1993). There are two ways such a discrepancy might be 
explained. The first is that the line has been smeared out 
greatly by Doppler motions. Another is that a large frac- 
tion of the fluorescent Ka line photons have been resonantly 
trapped and then destroyed by the Auger effect (sect. II of 
Ross, Weaver & McCray 1978; Ross & Fabian 1993, here- 
after RF). This process requires that the iron atoms being 
photoionized have between four and ten electrons; i.e. be in 
the form Fe xvn-Fe xxm. Comptonization is another mech- 
anism that can destroy the line, but the temperature and 
optical depth of the Comptonizing medium must be such 
that it destroys the line but not the edge; i.e. kT ~ 5 keV 
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and Tr < 1. If a corona above the accretion disc provides 
the illumination then it must have kT 3> 5 keV to obtain 
the observed spectral shape and thus if it destroys the line 
it will destroy the edge as well. The corona in the model 
for Cygnus X-l of Haardt et al. (1993) with kT « 150 keV 
and tt « 0.3 will only destroy the line if the inclination of 
the disc exceeds « 70 degrees so the effective optical depth 
along our line of sight exceeds unity. 

The smeared iron K edges of BHC are much deeper 
than those expected and seen in cold AGN discs (cf. fig. 9 
of Matt, Perola & Piro 1991). This would naturally occur 
in the reflection model if BHC discs were more highly ion- 
ized than AGN discs, since stripping of low atomic number 
elements greatly reduces the photoelectric absorption below 
7 keV thereby making the iron edge in the reflection spec- 
trum much more prominent. Smeared iron edges that start 
around 7.1 keV, the edge energy for Fe I, need not imply 
that the iron in the disc is in a low ionization state as we 
shall show below (cf. Marshall et al. 1993). 

In this paper we compute spectra of accretion discs 
around BHC using the Fokker-Planck optically thick radia- 
tive transfer code of RF. Our radiative transfer and ioniza- 
tion structure calculations are carried out for fully non-LTE 
disc plasma (electrons are taken to be in a local Maxwellian 
distribution due their large elastic scattering cross sections 
but the radiation field is not Planckian and the atoms do 
not follow the Saha-Boltzmann law) and incoherent Comp- 
ton scattering is included. We exhibit the predicted smeared 
edge profiles from accretion discs and discuss how they can 
be used to learn about disc properties. We also detail the 
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resonant Auger destruction process and point out its impli- 
cations. The results presented here, which highlight the edge 
and line profile in BHC, are complementary to and extend 
those presented in Matt, Fabian & Ross (1993) which con- 
centrate on the emission line profiles from irradiated discs 
in AGN. 



2 THE MODEL AND COMPUTED DISC 
SPECTRA 

The geometry we adopt is that of X-ray reflection and is sim- 
ilar to that assumed in much previous work (e.g. Ebisawa 
1991; Done et al. 1992; Marshall et al. 1993; Ueda et al. 
1994). Simple fits of the iron absorption edge suggest that 
any absorbing gas is only moderately ionized and therefore 
any transmission model involves large amounts of soft X- 
ray absorption. This is contrary to observation and so sim- 
ple transmission models can be ruled out. Partial-covering 
models, in which some of the underlying power-law contin- 
uum escapes directly to us, alleviate this problem somewhat 
and also provide acceptable fits (Marshall et al. 1993). 

In reflection models the primary power-law continuum 
is both directly visible to us and incident on plasma where 
it is absorbed, scattered and reprocessed to yield a further 
observed component, the reflection spectrum. The iron edge 
and any line would be associated with this component. The 
plasma is usually assumed to be the inner parts of the ac- 
cretion disc around the black hole. 

Our basic model is that of a geometrically thin, a- 
viscosity accretion disc (Shakura & Sunyaev 1973) around 
a 10 M Schwarzschild black hole. We consider the disc re- 
gion between 6 and 26 black hole gravitational radii (r g = 
Gm/c 2 ) and consider the radiative transfer in the outer 3- 
10 Thomson depths (this corresponds to a physical depth of 
~ 50 m). We adopt cosmic element abundances (Morrison 
& McCammon 1983) and energy generation rate propor- 
tional to gas density. The resulting equations for the hydro- 
gen number density (hh) and the disc flux (F s ), valid for 
a radiation-pressure-dominated disc, can be found in Ross, 
Fabian & Mineshige (1992). We take the flux distribution 
deep in the disc to be a blackbody with temperature T s so 
that B V (T S ) = F s but then solve for the full radiative transfer 
of the photons in the outer Thomson depths, including in- 
coherent Compton scattering. We take the viscosity param- 
eter a = 0.1 and the energy conversion efficiency r\ = 0.083 
(Shakura & Sunayev 1973). 

We consider disc illumination by an optically thin 
corona in which magnetic dissipation takes place, as is illus- 
trated in fig. 3 of Galeev, Rosner & Vaiana (1979). We as- 
sume Compton, bremsstrahlung, and synchrotron processes 
in the corona cause it to radiate a power-law spectrum with 
photon index T = 1.7 which isotropically strikes the disc 
from above and photoionizes its surface layers. We take the 
coronal illuminating flux at the surface of the disc, Fh, to 
equal F s (a reasonable assumption for lack of any alterna- 
tive since the illumination power is ultimately derived via 
accretion). Details of the treatment of the illumination are 
given in RF. 

Using the code of RF, we compute spectra from 10 irra- 
diated annuli between 6-26r g , chosen so that the flux from 
each annulus is roughly 10 per cent of the total flux con- 



sidered. The successful operation of the code requires the 
plasma to be highly ionized which is why we do not consider 
radii greater than 26r g or luminosities less than 0.1Z/Bdd- 
Line formation is calculated in an escape probability for- 
malism (Ross 1978), and we assume that fluorescence by 
FeXVII-FeXXIII is suppressed by resonant scattering (see 
Section 4). We add the annuli fluxes together taking into 
account Doppler blurring and gravitational redshift as de- 
scribed Chen, Halpern & Filippenko (1989). We have ig- 
nored the effects of light bending and our results should 
strictly be interpreted as time averaged since we do not take 
into account time dependence of the radiative transport, 
thermal transport, or illumination. A comparison of rapid 
X-ray variability (~ 1 ms; Rothschild et al. 1974), reflection 
region light crossing (~ 26r g /c ~ 1 ms), and plasma recom- 
bination/heating/cooling timescales (<C 1 ms) shows that at 
any instant disc plasma will be a 'quasi-thermal mess' state 
where electron temperature and ionization structure vary in 
a complicated way over the disc surface (note 'flickering' il- 
lumination need not remain fixed in one place) . We also add 
the outward flux from the corona assuming it is an isotropic 
radiator. 

Fig. 1 shows our results for discs at an inclination angle 
of 45 degrees. The total flux considered is roughly 77 per 
cent of the total flux from the disc plus corona system if 
the corona illumination is taken to extend to « 26r g . Just 
above 7 keV a smeared edge feature is formed exactly as 
reported by Ebisawa (1991). Iron Ka lines are conspicuous 
only in their absence; resonant Auger destruction destroys 
a large fraction of line photons and Doppler motions smear 
the rest from visibility. The low energy shape of the smeared 
edge is influenced by line photons which tend to fill it. The 
peaked structure around 0.7 keV in the 0.1Z/Edd case arises 
primarily from iron L lines (mainly from Fe xviii-xxi), but 
there is also a contribution from the K line of O vm. This 
structure may be compared with that shown in fig. 1 of Hess 
et al. (1994). 

The contributions from some of the annuli in the disc 
operating at O.lLBdd are shown in Fig. 2 (weighting by annu- 
lus area has been included). Several of the innermost zones 
produce a strong emission line from highly ionized iron, but 
this is blurred beyond recognition if i 20 degrees when 
Doppler and gravitational effects are included. The outer 
zones have little line emission due to resonant Auger de- 
struction. When i £ 20 degrees, the emission line from the 
innermost radii becomes apparent. 

One important facet of our model is that an edge fea- 
ture can be formed around 7.1 keV even from iron that is not 
cold, principally because of Doppler blurring, gravitational 
redshift and transverse Doppler shifts. The edge energy can 
vary between 7 and 8 keV depending primarily on the in- 
clination and disc ionization parameter (its energy increases 
with both increasing inclination and increasing ionization 
parameter). This is consistent with current observations of 
the smeared edge (cf. sect 6.3 of Ebisawa 1991), and ASCA 
will allow a more precise test of this prediction of our model 
and reflection in BHC more generally. 
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Figure 1: BHC Disc Reflection Spectra 
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Figure 1. Spectra from the region 6-26r g of illuminated BHC discs at a 45 degree inclination angle around a 10 Mq black hole. Note 
the prominent smeared edges and the lack of a visible iron Ka line. 



3 COMPARISON WITH OBSERVED 
SPECTRA 

We have converted our O.ll/Bdd model spectum at i = 45 
degrees into a fake X-ray spectrum assuming the response 
of one of the Gas Imaging Spectrometers, specifically GIS2, 
on ASCA (cf. Tanaka, Inoue & Holt 1994) and fitted the 
data over the 4-10 keV band with several models using the 
spectral fitting package xspec (Shafer et al. 1991). As ex- 
pected a simple power-law is a very poor fit. To illustrate 
this we have fitted a power-law over the 4-6.8 keV band and 
extrapolated it to 10 keV. The ratio of the (fake) data to 
this model, shown in Fig. 3, clearly shows the beginning of 
the large broad iron edge and the lack of a line (the small 
deviations from the power-law below 6 keV are numerical 
artifacts of the discreteness of the annuli and the blurring 
processes). The maximum fractional depth, of about 23 per 
cent, is approximately twice as large as that observed in 
data of Cyg X-l (e.g. Ebisawa 1991). A ten per cent in- 
crease in the direct continuum seen relative to the reflected 
spectrum, perhaps due to a different inclination angle than 
that assumed here, could easily reduce the predicted edge 
to agree with observations. 

Ebisawa (1991) and Ebisawa et al. (1994) have modelled 
the absorption feature empirically with a 'smeared edge' 
model (denoted 'smedge' in xspec) in which the iron ab- 
sorption cross-section, a, above the edge energy i3 B d g e has 
the form a[l — exp(— (E — EEd ge )/Ew)]- We confirm that 
this gives an excellent fit to the fake data from our mod- 
els. At 0.1Z/Bdd we find that a power-law with photon index 
1.45, £ B d g e = 6.93 keV, r max = 1.23 and E w = 3.72 keV 



yields a fit in which the residuals are all less than 1.2 per 
cent of the data. r max is the maximum optical depth in the 
smeared edge model. Ebisawa (1991) reports i?Ed g e ~ 7keV 
and Ew > 4keV from Ginga data of several BHC. At this 
stage we conclude therefore that for intermediate inclina- 
tion angles the blurred absorption feature produced by our 
model has a similar shape to that observed from BHC, a 
success for the reflection model. 

To assess the influence of resonant Auger destruction on 
the final iron Ka line observed, we have constructed further 
fake data in which the disc in the range 6-26r g is assumed to 
emit the spectrum appropriate for the zone centred at llr g , 
where the iron is principally FeXXV (Fig. 2) and Auger 
ionization cannot occur. A broad emission feature is then 
clearly seen at energies below the edge. Fitting these fake 
data with a model consisting of a power-law, a smeared edge 
and a Gaussian line yields an equivalent width for that line 
of 390 eV (the r.m.s. line width is 0.9 keV). This contrasts 
with the situation in our basic model at O.lLBdd where no 
line is required (see Fig. 3); any line (which is no so broad 
as to encompass the edge) has an equivalent width less than 
15 eV. This suggests that Auger destruction is an impor- 
tant reason why a detectable iron line is absent in our basic 
model. Doppler blurring at small radii is the other reason. 



4 RESONANT AUGER DESTRUCTION 

We now turn to an examination of the process of resonant 
Auger destruction of iron Ka photons and its effects on 
the iron line emitted by a highly photoionized accretion 
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Figure 2: Annuli spectra 
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Figure 2. Contributions to the iron Ka line from annuli centered 
at the radii shown for a disc operating at O.lLgdd- The predomi- 
nant iron ions at the surface of the disc as a function of radius are 
Fe xxv (8.6r g ), Fe xxv (11.0r g ), Fe xxv (14.0r g ), Fe xxi (18.0r g ) 
and Fe xix (24.0r g ). 
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Figure 3. The data to model ratio when a power-law model is fit 
to faked ASCA GIS2 data (made from our 0.1I/Bdd computation 
viewed at 45 degrees) in the 4-6.8 keV band and then extrapo- 
lated to 10 keV. The beginning of the large broad iron K edge at 
~ 7 keV and the lack of a line are apparent. 



disc. The detailed atomic physics of iron line formation is 
formidable as is illustrated by high resolution solar flare and 
tokamak spectra (e.g. Hill et al. 1979; Phillips et al. 1983; 
fig. 1 of Seely, Feldman & Safronova 1986; fig. 7 of Jacobs 
et al. 1989; fig. 5 of Beiersdorfer et al. 1993). We will of 
necessity restrict the scope of our treatment. 

We consider iron line formation in an accretion disc 
plasma with n e « 2 x 10 21 cm~ 3 and T e « 5 x 10 6 K. 
Many iron Ka photons produced by fluorescence of Fe xvn- 
Fexxn will be destroyed by the combination of resonant 
trapping and the Auger effect. The Ka fluorescence lines 



for each of these iron species result from transitions to lev- 
els of the ground electronic configuration of the next more 
highly ionized species. This next more highly ionized species 
coexists with the fluorescing ion in the photoionized plasma, 
and collisional excitation establishes a statistical population 
among the various levels of its ground electronic configura- 
tion. Therefore these lines are strong resonance lines, and 
the Ka photons are easily destroyed via reabsorption in the 
line followed by autoionization (Auger effect). 

To clarify the situation, consider the specific case of 
Ka fluorescence by carbon-like Fe xxi. A K-shell ionization 
of Fe xxi in its ground electronic configuration leaves the 
ls2s 2 2p 2 configuration of Fe xxn. Jacobs et al. (1989) have 
listed the eight most important radiative transitions from 
this excited configuration to the ground ls 2 2s 2 2p configu- 
ration, along with calculated transition rates and branch- 
ing ratios. Each transition has a spontaneous emission rate 
A n i ~ 10 14 s _1 . assuming a Doppler line profile, the aver- 
age opacity in one of these lines compared to the Thomson 
scattering opacity is 

Kline _ ntjne 2 /m e c)ft u 

KT ?l e (TTAVD 

where n; is the number density of ions in the lower level 
(either the 2 Pi/ 2 or 2 P 3 / 2 ground level of Fe xxn), n e is 
the number density of free electrons, fi u is the absorption 
oscillator strength, and Ai^d is the Doppler width. Adopting 
the abundances of Morrison & McCammon (1983), we get 

where n^e is the total number density of iron ions, A is the 
wavelength of the line, and g u and gi are the degeneracies 
of the upper and lower levels, respectively. Even if the lo- 
cal fraction, ni/n^e, of iron ions in the appropriate ground 
level is only a few per cent, the line opacity greatly exceeds 
the Thomson opacity. Since iron fluorescence occurs within 
a layer of Thomson depth tt ~ 1 (see George & Fabian 1991 
and references therein), resonant trapping of the line pho- 
tons is important. Following the original K-shell ionization 
or a resonant reabsorption, the probability of a radiative de- 
excitation (the 'fluorescence yield' V) varies from a low of 
Y ss 0.1 to a high of Y « 0.8, depending on the particu- 
lar excited level involved. (Jacobs & Rozsnyai 1986 found 
an average fluorescence yield Y = 0.48 for Fe xxi). Even 
for transitions with the highest fluorescence yields, only a 
few reabsorptions are needed to destroy most of the Ka line 
photons. A line photon is destroyed after diffusing only a 
fraction of a Thomson mean free path away from its origi- 
nal point of emission, and few line photons escape from the 
gas. 

Fluorescence by beryllium-like Fe xxm is inherently 
weak. Single-photon transitions would be expected to be 
unimportant for the ls2s 2 configuration following K-shell 
ionization. Although Chen et al. (1981) found that configu- 
ration interaction with ls2p 2 increases the fluorescence yield 
substantially, it is still only Y = 0.11. Resonant Auger de- 
struction can reduce the yield even further if 2s — > 2p colli- 
sional excitations in the dense plasma produce a large popu- 
lation of Fe xxiv ions in the ls 2 2p configuration (see Jacobs 
et al. 1989). 

The fact that Auger destruction predicts a substan- 
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tial weakening of Ka lines from intermediate iron ionization 
stages suggests that it may be part of the solution to what 
would otherwise be an unresolved problem: Fe xvm-Fe xxiv 
Ka lines (which result from photoionization of Fe xvn- 
Fe xxm) have never been observed from Thomson thick 
plasma in a cosmic X-ray source. (This is certainly the case 
for BHC, and we showed in Section 3 that Doppler blurring 
by itself has difficulty explaining the lack of iron lines in the 
reflection model for BHC spectra.) We encourage further re- 
search on the detailed treatment of iron Ka transitions in 
the context of Thomson thick radiative transfer in order to 
address this important issue. 

Another point worth noting is that bulk turbulent mo- 
tions in the direction perpendicular to the disc plane might 
reduce the line opacities enough so that Fe xvm-Fe xxiv line 
formation can take place. Bulk turbulent velocities larger 
than c times the ratio of the line width (~ 1 eV) to the 
line energy would be required (about 40 km s _1 ). Turbulent 
velocities of ~ ac SOU nd depend on the viscosity law chosen, 
and would be ~ 30 km s _1 if Cg 0und = P gas /p (cf. sect. 7 of 
Pringle 1981; we use a = 0.1 as per Sect. 2). Note in theory 
measurements of or upper limits on the Fe xvm-Fe xxiv 
lines from BHC might be used to learn about the viscosity 
law in the inner parts of their discs. 
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